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1 . 


INTRODUCTION AND SUMMARY 


Two measures of the relative performance of satellite based data 
collection and location systems are the precision with which platforms can be 
located (and/or their velocity estimated) and system capacity (average bit rate 
handled by the satellite) . This report investigates two concepts whereby im- 
in location precision and system capacity can be achieved relative 
to current systems . One of these concepts (TRILOC) can provide either or both 
Improvements in location and system capacity by utilizing three measurements 
acquired during a single transmission from a platform for location pxaposes. 

The other concept is to Increase system capacity by means of directive antennas 
on board the satellite for the purpose of reducing interference (either between 


platforms or from RFI) and/or reducing platform power. 
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TRILOC concept is evaluated in Section 2 , The location algorithm 
is developed ^as sum ingV ange -rate . radial-acceleration, and range differences 
are measured on-board the satellite during a transmission from a platform. Ar. 
error analysis of this algorithm is then derived from which platform location 
errors are related to the precision with which the, th ree measur ements are made. 
An example TRILOC processor is then postulated ;to^^^s^ m['a te |he threg^ measure- 
ment precisions.] 




Particular conclusions related to the TRILOC concept caVr be summ- 


arized as follows 


• If all measurements contribute equally to the precision 
of the location estimate, the improvement in location 
precision is basically due to a. greater number of measure- 
ments being available for noise suppression rather than 
an inherent advantage associated with combinations of 
different types of measurement. 
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• From ths analysis of th© ©xampl© TRILOC procossor, 
the errors in measuring radial acceleration cause much 
larger location errors than either measurements of 
range -rate or range difference. 

o The value of range rate and radial acceleration measure- 
ments decrease rapidly with increased satellite altitude 
(e.g. from a two to three hour period) while the value of 
range difference measurements does not change appreci- 
ably. 

An overall characteristic of a TRILOC system is the opportunity to 
improve location precision without sacrificing system capacity or alternatively, 
to Increase system capacity without degrading location precision. If system 
capacity is dominant, then the acquisition of three measurements of location 
parameters during a single transmission provides a nearly threefold Increase in 
system capacity compared to current data collection and location systems . On 
the other hand, if precise location is most important, then the additional loca-j 
tion parameter measurements acquired during each transmission can be utilized 
to Improve location accuracy without reducing system capacity of current 
systems . 

The potential advantages and disadvantages of both sweeping and 
fixed directive antennas on— board the satellite for random access systems are 
analyzed in Section 3 . These analyses indicate that there is no advantage 
(and probably only disadvantages) in having a sweeping antenna on-board the 
satellite compared to multiple directive antennas from either a performance or a 
system capacity viewpoint. This conclusion is based upon a necessary in- 
crease in mutual interference between platform transmissions as the gain of 
a sweeping antenna increases. However, directive antennas do offer advantages 
in system capacity improvement by decreasing interference , platform power 
reduction, and/or RFl suppression. 
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To achieve the advantages of the directive antennas, the beams must 
be of elliptical cross section and maintained in orientation relative to the 
satellite's direction of motion. In particular, the major axes of the beams 
must be parallel to the satellite's sub-track and be able to receive trans- 
missions from platforms from one edge of the satellite's visibility circle to 
the other. System capacity will be directly proportional to the number of 
these beams (from an interference standpoint) that are used to completely 
illiminate the visibility clrcle-l.e. the minor axes of the beams determine 
the number of beams required and also the improvement factor. 

To determine antenna size, the minor axis of the beams and 
operating frequency must be specified. To estimate required sizes for a 
low altitude satellite (Nimbus or TIROS), a synthetic aperature or array is 
assumed at an operating frequency of 400 MHz. This leads to the approxi- 
mate relationship that the maximum antenna dimension in meters la about .32 
times the number of beams. 
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2 . 


TRILOC CONCEPT 


2.1 Position Determination 

The bases for satellite referenced location systems Is the accxorate 
knowledge of satellite position and velocity throughout its orbit and the ability 
to measure range or range rate between the satellite and the platform to be 
located. For example, in the case of the IRI5 system, range is measured at 
two points during an overpass of the satellite. At each point, this measured 
range plus known satellite position establishes a sphere centered at the 
satellite upon which the platform is known to lie. Assuming the platform lies 
on the earth's surface, then the two range spheres and the earth's sui*face have 
two points of simultaneous intersection— one of which conres ponds to the geo- 
graphic coordinates of the platform . 


An entirely analogous location technique is that utilized by the RAMS 
system. Here, the measured quantity between satellite and platform is tange- 
rate which give rise to cones with apexes at the satellite (instead of spheres). 
The location is accomplished in the same manner, however, by determining the 
two points of simultaneous intersection of two cones and the earth's surface. 


A third example is the TRANSIT navigation system that measures the differences 

in range between satellite and platform to estimate location. 

The TRILOC concept is to measure not Jus ^ rang^ (or ^ra^e->dif f e r- 
lence ) .n^/ranoe-rat^hnt to measure both of these quan^'esandCqg ^al a cceleratl^ 
dur ing any sin gle t ransmission from a platform. For example, if range and range- 
rate were measured during a single transmission from a platform, then its loca- 
tion could be estimated by determining the two points of simultaneous intersection 
of a sphere (range), a cone (range-rate) and the earth's surface (assumed 
altitude of the platform). Based on this alone then, addition of a range measure- 
ment during each transmission of a RAMS platform would effectively halve the 
number of transmissions required for location purposes. 
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As mentioned, radial acceleration Is also capable of being measured 
and is as equally useful as range or range-rate for location purposes. In this 
case , the geometric surface upon which the platform is known to lie can be 
shown to be shaped like a doughnut with no center hole and whose center is 
at the satellite’s position when the measurement is made. 

Conceptually then, as many as three measurements could be obtained 
from a single platform transmission for purposes of location. Whether this is 
possible, particularly, in the case of the range measurement, depends upon the 
communication link between the satellite and platform— i.e. , range measurement 
requires a two way link, as in the IRLS system. On the other hand, measurement 
of range difference between two transmission points can be accomplished with 
one way communication links. 

2 .2 Position and Velocity Determination 

The above discussions indicate platform position can be determined 
with at most half the number of platform transmissions of current location sys- 
tems by utilizing the TRILOC concept. The more general need, however, 
{particularly when Increased precision of location is required) is the ability to 
estimate platform velocity as well as location, Fxirthermore , precise location 
requires correction for center frequency drift and perhaps drift rate during a 
satellite overpass. Instead of two measurements being required, then, precise 
location probably requires four to six measurements during an overpass (if not 
more to permit noise suppression filtering) . Because three measurements are 
acquired for each transmission received In a TRILOC system, the necessary four 
to six measurements can be obtained with as few as two transmissions per over- 
pass . Therefore , implementation of the TRILOC concept can result in a factor 
of three reduction In number of transmissions required per overpass . 








2.3 LocattonA’elocity Estimation 

^ ^ ^ While a geometric Interpretation of satellite-referenced location 

^,iisc,ilitates description of the TRILOC concept, the actual method of deriving 1 
I location and velocity estimates is accomplished in a different manner. /ihis 
will be described first by limiting discussion to a system Ignoring platform 
velocity and relying solely on range measurements to estimate location. 
Secondly, the description will be extended to a TRILOC system ignoring platform 
velocity. These discussions will then be used as the basis for describing 
a TRILOC system estimating position and velocity while correcting for 
long term drift in platform transmitted frequency. 


2.3.1 Range Measurement Location 

If two measurements of range between the satellite and platform are 

• 

obtained, the above discussions Indicate the location of the platform to be one 
of two points of simultaneous Intersection between two range-spheres and the 
earth's surface. Another way to establish the platforms location Is to compute, 
a priori, for each possible location of the platform the two corresponding values 
of range that would have been measured if, in fact,, the platform were at that 
location. When the actual pair of range measurements Is acquired, the location 
of the platform can then be determined by matching the measured pair to a com- 
puted pair. Note, except for a usually resolveable ambiguity, only one of the 
computed pair will match the measured pair. 

Actually the necessity to pre-compute all possible pairs of range 
measurements isn't necessary if the following procedure is followed. First, 
establish an estimate of the geographic coordinates of the platform (X, V) . 

This estimate can be as crude as requiring the position to be within view of the 
satellite when both of the range measurements are made. Secondly, for this 
estimated position and the position of the satellite when the actual range mea- 
surement were made, compute the pair of measured ranges that would have 
occurred if the estimated location were correct. These will of course differ from 
the measured ranges— unless the guessed location Is coirect. 
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The difference between the measured and computed ranges are then 
used to correct the X,Y estimate of platform location to more closely correspond 
to the actual platform location. This is accomplished by equating the difference 
in each case to the sum of two terms • One of these terms is that portion of the 
(measured~computed)range difference occurring because the guess of the X 
coordinate is in error and the other term is the Increment caused by the error 
in the estimated Y coordinate of the platform. Symbolically, this may be 
written as 




where fix and 6y are the errors in the X and Y estimates of platform location 


is the difference in(measured minus computed) range that 
will occur at the first measurement point if the position 
estimate were In error by AX 

(ARj/AY), (AR3/AX) — — analogous to (ARi/AX) 

fs that portion of^Rj measured - Ricomputed) caused by a 6x error 
in estimating the X coordinate of platform location. 

However, if these two equations are solved simultaneously for 6X and 
fiY and these values are used to correct the original X,Y estimate of platform loca- 
tion, then the result will be a better estimate if not the actual location of the platform. 

An Important extension of this type of solution process concerns those 
conditions wherein the number of range measurements is greater than two— 
i.e. , the number of measurements is greater than the number of unknowns. If ’ 
there are no errors present in the range measurements or satellite position, 
such redundant measurements cannot of course serve any purpose. However, 
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with errors present / redundant measurements can be utilized to suppress the. 
effects of the errors provided they are random with negligible bias • This is 
accomplished by forming a sum of the squares of the differences between 
measured and computed (i.e. , based on estimated location) ranges at each 

♦ ^ 

measurement point and then determining the location of the platform that 
minimizes this sum. For the case where the errors are normally distributed, 
this process can be shown to decrease the error in locating the platform in 
proportion to the reciprocal of the square root of the number of measurements . 

2.3.2 TRILOC Location 

Analogous to the above, assume a TRILOC system acquires range and 
range-rate measurements between the satellite and platfprm. For the case wherp 
two measurements are acquired, the location process described above for a two 
range-measurement process is equally applicable. However, the extension 
to the redundant measurement solution is not quite as straight forward. 

In the case with redundant measurements of the same typef— I.e., 
all range or all range-rate) minimizing the sum of the squares of measured 
minus computed difference is conceptually, valid. However, for the case 
where different measurements are Involyedf—e.g . , range and range-rate — ) 
minimizing the algebraic sum of the squares of differences is not valid without 
assigning relative weights to the different measurements, let alone adjusting 
units to make the sum meaningful. 

One means of accomplishing this weighting of the different measurements 
Is to divide each difference by the standard deviation of the error for that measure- 
ment. This will, in effect, magnify those differences that correspond to low 
error measurements relative. to those of high error— additionally , the sum of the 
squares of these quotients becomes the sum of non-dimensional terms— i.e. , 
compatible units. This enhancement of the more precise measurement differ- 
ences serves to drive the location solution to that indicated by these measurements 
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Note this also means that when the standard deviation of one type of measure- 
ment error Is large relative to others, then this measurement is effectively 
ignored in the solution process . This can be seen in the data presented in 
Figure 1. 

Figure 1 is a presentation of the errors in measuring range differ- 
ence between successive transmissions from a platform, range rate# and 
acceleration normalized by the resulting mean location error-i,e. , the average 
magnitude of the location error. A basic assumption utilized to develop these 
data is that the measurement errors are normally distributed with zero mean. 
Furthermore, there are no errors present except the measurement errors. 

These data correspond to a satellite orbit similar to the Nimbus spacecraft. 
Additionally, the measurements are acquired over two successive overpasses of 
the satellite wherein the transmissions are separated by five minutes. The 
geometry of the overpass Is such that the sub-tracks are evenly spaced to each 
side of the platform and under these conditions, a total of 3 transmissions are 
received on each overpass. Note, this means a total of 3 measurements of 
range-rate and acceleration are received per overpass, but only 2 measurements 
of range difference are obtalned,for a total of 8 measurements per overpass . 

The first point of Interest for these data are the asymptotes noted In 

Figure 1 . These boundaries correspond to the condition wherein the contribution 

of two of the three measurements to the location estimate are Ignored because 

the standard deviations for these measurements are large relative to the third. 

Note, the asymptote for large errors In range difference and range-rate cannot 

be shown In Figure 1 , For example, if the standard deviations of range-rate 

and acceleration are large, then the standard deviation of range difference 

measurements establishes location error— i.e., mean location =a^j^/.S6. 

Similarly, if range difference and acceleration errors are large, then the loca- 

—3 

tion error boundary is <7^/3.48*10 

A more important result obtainable from Figure 1 , however, is the re- 
duction in mean location error when both range difference and range-rate 

measurements are used — assuming acceleration is not used. 
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Standard Deviation of Range Difference Measurements' 
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To demonstrate the advantage, assume cr^ is equal to «56 kilometers 

and is equal to 3.48 meters/second. For the asymptotes discussed above, 

R 

the mean location error, u. _ _ , for systems using one or the other of range-rate 

LOG 

or range difference measurements is found from; 


• Range -rate*-*the asymptote is described by 

a*/u = ,00348 per second. With d' equal to 
R LOG R 

- ,3.48 meters /s econd , is then 3.48/.00348 or 

1 kilometer. 


Range -difference —the asymptote is described by 

a 


'Ar 
^LOC 


equal to .56 kilometers. 


p is then .56/. 56 or 1 kilometer. 


To determine the mean location error when both measurements are used simul- * 
taneously, manipulation of the data provided in Figure 1 is necessary. 


With OP and cr • fixed at .56 kilometers and 3.48 meters/second and 
AR R 

with large (or Ignored) , the point on Figure 1 that establishes mean loca- 
tion error is the intersection of two lines . One of these lines is the = oo 
line of Figure 1 . The other line is one which must be drawn and is that line 
along which the ratio between and is a constant— In this case, 

O' = .56/. 00348 160 seconds . This line is sketched on Figure 1 and the 

R 

intersection between this line and the op = oo line is noted by the symbol ® . 
At this point, the ordinate and abcissa are; 


ordinate - 

abcissa — = 4.6 per second 

With either the ,56 value of or the 3 .48 value of the mean location 
error is then determined to be about .75 kilometers. 
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In the asymptotic case using range~rate measurements « a total of 
six measurements acquired during two successive overpasses are used. By 
adding the range difference measurements# of which there are four, there is a 
total of ten measurements being utilized. If the assumption is made that the 
reciprocal of the square root of the number of measurements determines location 
error, then the location error would be one kilometer divided by the square root 
of 10/6 or .77 kilometers. Because this is nearly equal to the .75 kilometer 
error derived from Figure 1 , this indicates that the number of measurements is 
the major cause of error reduction Instead of an inherent advantage associated 
with different types of measurements , However, it should be noted that if the 
number of transmissions from a platform are fixed, then the location precision 
is improved because of the multiple measurements acquired during each 
transmission. 

% 

0 

As another example , the location error resulting when all three 

measurements are used can be determined in a manner analogous to the above . 

The assumption is again made that each measurement has a standard deviation 

of error such that if it alone were used for location, the mean location error 

would be 1 kilometer. The standard deviation of the acceleration error Is 

therefore as noted in Figure 1“— .85x10 /sec. The standard deviations for 

range difference and range-rate are as before ,56 and 3.48 x lo"^/sec 

respectively. With these values, the location error is found to be about .6 

kilometers. This is established by the intersection of the or .^/cr . = 160 sec 

/SR R 

line and a curve (not shown) of constant cr.. /tr. = .$5 x 10 “®/.00348 - 

R R 

.0025 sec. This corresponds to the use of 16 measurements-— six range-rate, 
six acceleration, and four range difference. Note again that 1// 16/6 also 
gives an estimated error of about -6 kilometers. 

2.3.3 TRILOC Performance 

Extension of the previous discussions to the case wherein platform 
velocity and mean frequency of the platfcrm transmitter during an overpass are 
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also estimated is not significantly different.* Instead of the relationships 
given in Section 2.3.1 for X, Y determination only, each (measured -computed) 
difference is now equated to the sum of five terms corresponding to two velocity 
coordinates (e.g. , east and north \mder the assumption of negligible or known 
vertical velocity) , mean transmitted frequency and again the two position 
coordinates. Additionally, a parameter held constant In previous performance 
descriptions is satellite altitude. Because of the potential variations in satel- 
lite altitudes , the effects of this parameter should also be evaluated. There- 
fore, discussions and data provided below are concerned with TRILOC 
perfcarmance in two ways; 

• The variation of platform location error with the 
standard deviations of range difference, range- 
rate, and acceleration measurements—the units 

A 

in this case being anticipated units of the actual 
measurements to be made, namely micro-seconds for 
range difference , Hertz for range-rate , and Hertz/ 
second for acceleration 

• The variation of location error with satellite altitude, 

r 

holding the number of platform transmissions con- 
stant by Increasing the interval between transmissions 
for higher altitude orbits . 

To evaluate the effects of satellite altitude on the location errors of 
a TRILOC system, three altitudes are considered for sun-synchronous orbits— 

108 minute, 2 hours, and 3 hours. For these orbits, the number of trans- 

t 

missions during which range -difference, range-rate and acceleration are measured 

* Section 3.4 provides the mathematical derivation in detail, while Section 
3.5 analytically describes the means by which location error is computed 
as a function of measurement errors. 
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Is held fixed at four per overpass with two overpasses . This means there 
are a total of four acceleration# four range rate and three range difference 
measurements per pass or a total of 22 measurements over the two successive 
overpasses . The overpass geometry for these three orbits are sketched In 
Figures 2 # 3 # and 4 « 

In Figure 2 , the overpass geometry for the 108 minute orbit is shown. 

The platform to be located is taken to be at 0® latitude and 0° longitude and , 
for the computations of the location error, has zero velocity. As noted, the in- 
clination of the orbit is 100^. The successive passes are seen to be symmetri- 
cally placed east and west of the platform position and the circles denote the 
satellite sub-point at the eight transmissions. The transmissions are separated 
by four minutes. Figures 3 and 4 provide similar information for the geometries 
of the two and three hour orbits. For the conditions of these figures, little if any 
geometric dilution of precision should occur in deriving location of the platform. 

Figures 5, 6, and 7 present the variation of location error for a TRILOC . 
system operating with the three different orbit altitudes and geometries des- 
cribed in Figures 2,3, and 4. As previously mentioned, these data correspond 
to estimation of two position coordinates , two velocity coordinates and the 
mean transmission frequency of the platform. The presentation Is entirely the 
same as that described in Figure 1 when only platform position coordinates 
were estimated from the measiirements (except for the introduction of the units . - 
noted) . 

As a first point of interpretation, the improved location accuracy of 
a TRILOC system can be derived from these Figures at the three different 
altitudes. This is accomplished by assuming the standard deviations of the 
measurement errors are equal to those for the single-type-measurement system, 
that has the performance of a one kilometer mean location error. For example, • 
in the case of the 108 minute orbit. Figure 5, these standard deviation asymp- 
totes are .89 hertz for frequency (range-rate) measurements, and .0116 Hertz/ 
second for frequency rate (acceleration) measurements. In this case, a TRILOC 
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Period » 120 Min (Altitude • 1680 Km) 


Degrees 

Latitude 



FIGURE 3 . OVERPASS GEOMETRY 
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Pe iod « 180 Min (Altitude - 4182 Km) 


Degrees 

Latitude 



FIGURE 4* OVERPASS GEOMETRY 
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TRILOC Location Errors 


# Solving for 

Position (Two Coordinates) 

Velocity (Two Coordinates) 

Transmission Frequency on Each Overpass 

• Overpass Geometry of Figure 2 
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Standard Deviation of Frequency Measurements 
Mean Location Error 
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FIGURE 5 
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Mandat'd Deviation of Time Difference Measurement I 
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TRILOC Location Error* 


• Solving for 

Position (Two Coordinates) 
Velocity (Two Coordinates) 
Frequency on each Overpass 

• Overpass Geometry of Figure 3 
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Standard Deviation of Freq uency Measurement 
Mean Location Error 
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FIGURE 6 
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TRJLOC Location Errors 


• Solving for 

Position (Two Coordinates) 
Velocity (Two Coordinates) 
Frequency on each Overpass 

• Overpass Geometry on Figure 4 



FIGURE 7 
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system would yield a location error of approximately .53 Km which again Is 
not significantly different than one kilometer divided by /22/8‘. 

The variation of location error with altitude may be determined by 
comparing TRILOC systems with equal measurement errors. For example, 
if the errors in the above examples are used to enter Figures 6 and 7 , then 
the location errors are found to be .66 Km for the two hour orbit and .81 Km 
for the three hour orbit. These should then be compared to the .53 Km loca- 
tion error for the 108 minute orbit. 

An interesting side light of Figures S through 7 concerns the varia- 
tions of the asymptotes — l.e, , the variation of location errors for 'those 
systems utilizing only one type of measurement. In the case of both range- 
rate and acceleration, increasing satellite altitude degrades location per- 
formance . For a range rate system , the required standard deviation of error 
for one Kilometer location error is 2,1 Hertz at 108 minutes and ,67 Hertz at 
3 hours— or, for given measurement error an increase by a factor of more 
than three in location error. For an acceleration only system, this facior 
of increase is about 14 . On the other hand, an increase in satellite al- 
titude fora range -difference-only system results in a decrease in location 
error. This discrepancy in location error trends with altitude for the single 
measurement systems can be qualitatively justified by the following argument. 

As the satellite altitude increases, the total variation of both range- 
rate and acceleration during an overpass will decrease. Furthermore, the 
variation in acceleration will decrease more rapidly with Altitude than range- 
rate. Therefore, holding. measurement errors constant result in relative pre- 
cision of measurement to the variation of the quantity being measured to de- 
crea se — i . e . , larger location errors . 
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This is not tnie for range difference and in fact, the reverse is 
true. In this case the total variation of range differences during an overpass 
increase in going from a 108 minute orbit to a two and three hour orbit. 
Therefore, the location error for given measurement error should decrease. 

A more important aspect of these variations to a TRILOC system is 
that the higher the satellite altitude, the less advantageous a TRILOC system is. 
The trends suggest that at some altitude, depending upon measurement 
errors that are achievable, the location error will be solely determined by 
the range difference errors-— l.e. , operation will be limited to the asymptote 
where the standard deviations of range-rate and acceleration are very large 
relative to the precision of range difference measurements . 

2.3.4 TRILOC System 

While the previous discussions indicate the advantages inherent In 
the multiple measurements of the TRILOC concept, estimates of attainable 
' location precision can only be made by postulating actual mechanisms for 
measuring the three parameters. To obtain these estimates, a system is 
suggested in the following discussions regarding the processer on-board 
a satellite capable of making the requisite measurements. From this, estimates 
of actual measurement errors are determined and the data of Figure 6 is used 
to evaluate location precision as a function of system capacity— l.e. , number 
of platforms. 

The signal spectrum of this example TRILOC system is a carrier with 
20 kHz sub-carriers that are modulated with data. The three parameters 
necessary for TRILOC location are measured during a transmission Interval as 
follows: 

• Received frequency and frequency rate are 

computed from a least square error regression 
of measured time Intervals between a successive 
number of positive zero crossings of the carrier 
while it is tracked by a narrow band phase lock 
loop 
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• The range difference measurement is performed by 
measuring the phase of the 20 kHz sub-carrier 
relative to a standard on-board the satellite, and 
comparing this phase to that measured during a suc- 
cessive transmission-phase ambiguities are not a problem 
In this case because the major portion of any change 
in range is due strictly to known satellite motion. 


To obtain estimates of the errors that can be anticipated in measuring 
these parameters and the resulting estimate in platform location precision, 
several assumptions will be made. 

e The significant system errors are those associated 
with the measurements due to noise 

♦ 

• The standards on-board the satellite do not con- 
tribute significant error to the measurements' 

• The overall signal power to noise density ratio 

Is 40 dB-Hz, equally distributed between the 

/■ 

carrier and the coherently demodulated subcarrier 

• The carrier Is tracked with a 20 Hz loop and the 
subcarrier is measured through a 400 Hz filter 

• The measurement period for the carrier Is ,5 
seconds and for the subcarrier is 5 milliseconds. 

With these assumptions, the signal to noise ratio for the frequency 
and frequency rate measurements will be 24 dB while the signal to noise ratio 
for the subcarrier measurement will be 11 dB.’ The phase noise in these two 
channels can then be determined from 


For carrier channel: phase noise variance 


2 S/N 


For subcarrier channel: phase noise variance = 


rs7N 


002 rad^ 

® .04 rad^ 
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utilizing the regression analysis provided in Appendix F , the standard 
deviations of the errors in estimating received frequency and its time rate 
of change are 

Standard deviation of frequency erron , .0022 Hz 

Standard deviation of frequency rate error: .0117 Hz/sec 

Where the assumption is made that the nominal 
frequency is 4 kHz and the period of every cycle 
is measured and used in the regression analysis 
to estimate frequency and frequency rate. 

Similarly, if the relative phases of 100 positive zero crossings of 
the 20 kHz subcarrier (5 milliseconds measurement period in order to be able 
to demodulate the data on the subcarrier during the majority of the transmission 
period) are measured and utilized to estimate the relative phase, then the 
standard deviation of the side tone phase measurement will be -/.04/100 or 
.02 radians. This translates into 

Range difference error =*/2 x range error = -fz x x wavelength 

For a 20 kHz subcarrier, the wave length Is approximately 15 km so that the 
standard deviations of the range difference error is approximately 70 meters . 

For purposes of utilizing the data of Figure 6, this is equivalent to .23 micro- 
seconds error (i.e. , the case where measurement of range difference was 
based upon signal transit time). 

These three errors .002 Hz, .01 Hz/sec and .23 microseconds—" 
coupled with the data presented in Figure 6 indicates the .002 Hz standard 
deviation in measuring received frequency will dominate in the location 
algorithm followed by the range difference measurement and lastly the frequency 
rate measurement. Therefore, Implementation of a TRUCK; system versus 
systems wherein only one parameter is used for location will depend upon the 
relative values assigned to precise location of platforms and the number of 
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simultaneous users in a system. Another equally Important consideration Is 
the ability to service a mixture of user platforms wherein location precision is 
Important for some and not for others . 

For those users requiring precise location, acquisition of a relatively 
large number of measxirements during an overpass is Important. If system 
capacity is limited by interference between platform transmissions, then this 
increased location precision necessarily results in a reduction in system 
capacity. This is particularly true for systems wherein only one location 
parameter is measured during each transmission. However, In a TRILOC sys- 
tem the degradation in system capacity is much slower. Compared to the single 
measurement system, a threefold increase in number of measurements is avail- 
able in the TRILOC system for the same number of transmissions . Furthermore, 
if the example processor previously described is Implemented, then the precision 
with which platforms can be located should be considerably improved compared 
to current systems typified by IRLS and RAMS. 

2.4 TRILOC Location Technique - Analytical 

The general problem of locating a platform may be formulated as 
follows: We want to compute a trajectory parameter vector x, whose elements 
may Include platform position, velocity, and acceleration, plus parameters 
relating to transmission frequency. Parameters relating to transmission fre- 
quency instabilities may include such items as the center frequency drift and 
the center frequency drift rate on each satellite overpass. We take a set of 

m observations , which will be denoted by the column vector z . In a Doppler 

th -* 

location system, for example, the j component of z would be the (observed) 
received frequency at tlmetj.. Finally, we assume the existence of an analy- 
tical model relating the trajectory parameter vector to the measurement vector, 
i.e., we assume that there is a known function F’such that F (x, t^) = Zj, the 


I 
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component of z . If the dimension (n) of x is equal to the dimension (m) 
of 2 , we seek the vector x for which 

F (x, tj) = zj U) 

for j = 1 , n . 

In practice, however, the dimension of z often exceeds the dimension 
of X, that is , there are more equations than unknowns , When this occurs , there 
usually is no vector x for which Equation (1) holds for all J because there are 
always some errors in the observations, and the analytical model represented 
by the function F may only approximate the physical system. When the number 
of equations exceeds the number of unknowns, the method of least squares 
may be used to compute x. 


The least squares technique consists of minimizinq the following 
quadraUc form; a, 






In equation (2) , F (x) is a column vector with components Fj “ F(x,tj), and 
is a weighting matrix, which is the covariance matrix of the measure- 
ment errors, A necessary condition for the quadratic form (2) to be a minimum 

- Ffx)J = 3 


where D is the differentiation operator 





The equations (3) are nonlinear and may be transcendental, making 
an analytical solution difficult or impossible. However, an iterative tech- 
mique (Newton-Raphson’s method) may be used to solve the system of Equa- 
tions (3) . If we have a sufficiently good estimate x^^ of the trajectory para- 
meter vector, we can expand F(x) into a Taylor series in a neighborhood of 

X. t 

F^X) “ l> ( F (X^)) f 

where we have retained only linear terms. 
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Equation (3) then become 




To simplify the notation, let A = D(fGco^) • A is then a matrix with 


elements 

= 

Equation (4) then becomes 




sx: 


or 


/p' iy~‘ [ z, F (Xc j ^ ~~ ^ 

fi']ix-xo) - [ (^'F(xc)) (^) 


In equation (5) , x - is the error in the estimate of the trajectory 

parameter vector, i.e, , the difference between the initial estimate and the 

true value of the vector, provided is sufficiently close to x. Similarly, 

2 “ F(x ) is the error in the measurement vector, i.e, , the difference between' 
o 

the actual measurement vector (z) and the value of the measurement vector 
computed from the initial estimate (>^) of the trajectory parameter vector. 

The solution may be iterated by beginning with the initial estimate 
Xq and computing successive improved estimates of x from 

x/t, = x: C ^ • F(Xi)) 


The solution is iterated until successive values of x remain 

\ V 

substantially unchanged. 


2.4.1 Selection of the Trajectory Parameter Set 

If three platform position components and three platform velocity com- 
ponents are to be computed at each time (t^) of transmission, then each trans- 
mission would have to supply at least six independent observables from which 
position/velocity can be computed. In most systems, only one observable is 
provided per transmission, and with the TRILOC system at most three are pro- 
vided, Therefore, a position/velocity estimate at time t^ cannot be obtained 


2-24 



independently of a similar estimate at some different time tj. Instead, a model 
of platform motion must be assumed which defines the relationship between 
platform kinematic parameters at different times* Formally, the model of plat- 
form kinematics must have the following form: 

Ag M ~ ^ ^ \ ^ •) 

These equations are interpreted to mean that if the kinematic state 
of the platform is known at time t^, then its state at any time t may be com- 
puted from the functions R, V, and A. In practice, the velocity and accelera- 
tion may be assumed to be zero, i.e« , the platform is assumed to be stationary. 
Alternatively, the platform may have a finite velocity but may always remain 
at a constant (known) altitude. It is this second case that is usually of interest 

A platform at constant altitude may be a stationary platform on the 
earth's surface, a buoy, or a constant altitude meteorological balloon. If the 
platform's altitude is constant, then one of its position components is a func- 
tion of the other two, and it has no vertical velocity component. If the plat- 
form follows a great circle route, then its acceleration is a function of its 
position and its velocity. Thus if the platform moves in a great circle route 
at constant altitude, the kinematic subset of the trajectory parameter vector 
contains only two components of position and two components of velocity. 
Specifically, these components are platfom position/velocity components at 
time tQ. These components are then used to compute the remaining position 
component, and the entire acceleration vector. Platform position/velocity/ 
acceleration at any other time t Is computed from tiie position/velocity/accel- 
eration at time t^ using Equations (6) . 
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In addition to initial position/velocity estimatss, the trajectory 
parameter vector may also contain -components identifying the platform trans- 
mission frequency instabilities . Platform location is computed by the informa- 
tion provided by the frequency received at the satellite. Location estimates 
can be no better than the Information on which they are based. Therefore, 
unless the characteristics of the transmitted signals are known or can be com- 
puted to a reasonable degree of precision, the estimates of position and velocity 
will be inadequate. The instabilities can be computed by assuming an unknown 
center frequency drift and, possibly, the center frequency drift rate. These 
unknowns are then computed just as the (unknown) position and velocity are 
computed . The center frequency drift and drift rate can be computed for each 
satellite overpass . * 

2 .4.2 The Observables and Their Relation to the Trajectory 
Parameter Vector 

i - 

* fl 

The TRILOC concept Includes the measurement three quantities. These 

are; 

1) the received frequency of a signal 

2) the time rate of change of the received frequency 

3) the difference between the time of reception of 
two signals (or successive measurements of side 
tone phase) . 

These observables are directly related to the relative kinematics between the 
satellite and the platform. 

The first observable is the received frequency, If the transmitted 
frequency is f^, then the Doppler frequency is 

4 = 4-4 
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The Doppler frequency is proportional to the time rate of change of the 
range between the platform and the satellite (V^)i 

e/ 






or. 


- h(>- 


o) 

In Equation (7) , c is the speed of light, and f^ is the actual transmitted fre- 
quency, which is computed from , 

where f^^ is the nominal oscillator center frequency, f^g is the center frequency 
bias at time t^, and f^ is the center frequency drift rate. is computed from 




it 




where R and V are the satellite position and velocity vectors at time t, and 

_» _S s ^ ^ 

Rg and Vg are the platform position and velocity vectors at time t. Rg and Vg 

are computed from the initial estimates of platform position and velocity and 

Equations (6 j. Equation (7), therefore, expresses the observable ^ as a function 

of the trajectory parameter vector, i.e. , all or any part of ttie elements of 

Rr / V„ , Ar , fr and f« Each Doppler measurement produces a relationship 
°o °o “o '■B '■ 

of the form of Equation (7) . 

The second observable is the time rate of change of the received 
frequency. Its relationship to the trajectory parameter vector is found by taking 
the time derivative of Equation (7): . 

h- o-fo-4-^ ■ (1) 

• • 

Again , ^t " ^to ^tB ^ ^t (t-to) , and c is the speed of light. A^ is the time 
rate of change of the relative range rate between the satellite and the plat- 
form, i.e. , ^ 

j ^ /fi “ ) •( — -f I ^ 

** ~ 1 ^ ^ , -A ^ 




Jt 
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Ag and Ag are the acceleration vectors of the satellite and the platform . 
Equation (8) then expresses the observable fg in terms of the components 
of the trajectory parameter vector. There is one equation of the form (8) 
for each TRILOC transmission. 


The third observable is the difference in time*between the recep- 
tion of transmissions( At) , The time required for a signal to travel from 
the platform to the satellite is proportional to the relative range (Rn) between 
the platform and the satellite, where the constant of proportionality is assumed to 
be the speed of light. If the time (A t^). between transmissions from the platform 
is known, then the difference in relative range ( ARj^) at two transmissions 
is proportional to the difference in time between the reception of the two 
signals, less the difference in time between the transmissions: 

tise of successive measurements of side tone phase achieves the same 
purpose with an entirely analogous derivation. 


A possible implementation of this measurement might be a platform that counts 
zero crossings of the signal provided by its oscillator and is programmed to send 
a transmission when the number of zero crossings reaches a predetermined num- ' 
ber Np. If the oscillator frequency is f^ then the time between transmission is 
Nc/f^. We then have 


or. 


AR{^ is computed from 




r?; 


»■ I ■“ I 

I 

-A ^ 

where Rg^ and RBi are the satellite and platform position vectors at time t^. 
Equation (9) then relates the observable At to the elements of the trajectory 
parameter vector. The number of equations of the form of Equation (9) on an 
overpass is one less than the number of transmissions on that overpass. This 
is because for the first transmission on an overpass there is no previous trans- 


mission from which At can be computed. 


2-28 


2.5 


TRILOC Position Error - Analytical 

In Section 2.4 the weighted least squares solution to the TRILOC 
location problem was derived as 

sx ~ CO 

For convenience, we will write x instead of <5x and & instead of 6 a. it is 
important to remember, hov/ever, that x and z are now the errors in the estimate 
of the trajectory parameter vector and in the measurement vector, not the 
vectors themselves. If the measurement errors are sufficiently small. Equa- 
tion (1) expresses the errors in the trajectory parameter vector as a function 
of the measurement errors . If f is the covariance matrix of the measure- 
ment error vector, then 

cov (K) ^ C [(/T'f-'/i) 

The joint probability density function for the vector x is 

j. 

4 tx) - -^—ji exP x'^Px ) 

*- 

where P = and n is the dimension of x. 


We want to find the expected value of the location error, i.e. , 
f 2 2 ~ ) 

E ^ (xi + X2 where we have re-arranged the elements of x (if necessary) 

so that XI and X 2 are the horizontal components of platform position. Then 
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The Integration is made easier by a transfonnation to hyper spherical 
coordinates, (defined in Appendix C) ; 

n~t 

!= CCSO, IT 
fil w 7T 

L 

^n-^i 

X/ = r 'Si' C t ^ c ^ n) 

7" A ^ 17 

{xU n - r CC4 

» 

k'^rx ^ r' ?*^p? 


The determinant / Jn/ of the transformation from (r, 0l , On-1^ coordinates 

to (xi , . . . , Xj^) coordinates is (see Appendix C) 

3 1 X t j » * * ) Vw 3 


I I / I y nj I /j-/ j[2~ ^ 

'■‘'•I “ * '' 7T 


9 




We then have 
£ 


'T T/x 


«o 
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Using the results from Appendices D and E, the number of integrations can be 
further reduced , so that 


£■ I Xx)"-] 



The symbols a 2 .6 2 » V * ^2 » ^2 defined in Appendix E. 

A closed form solution to the integral was not found, but its value 
was numerically computed to obtain the results of Figures 5 through 7 
previously discussed in Section 2* 
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3.0 


DIRECTIVE ANTENNAS 


' A directive antenna on-board the satellite may have distinct advantages 
for data collection and location systems relying on random access. Two such 
advantages would be the suppression of R FT and/or the ability to reduce radiated 
power from the platforms. Thirdly, becaus."- these systems are limited in 
capacity by Interference between platform transmissions, the use of a directive 
antenna may serve to alleviate this restriction also. 

The analyses that are required to determine the advantages and con- 
sequences of the use of a directive antenna (s) can be divided into three categories 

Sweeping Directive Antennas 

• Given that a specified number of transmissions 
from each platform is required during an overpass 
of the satellite , how are sweep rate and gain of 

? 

• « 

the antenna related to the.duty cycle of the platforms ? 

e Given the relationship between duty cycle and 

antenna gain, how does increased antenna gain 

effect the probability of mutual interference 

/ 

between platform transmissions? 

FTxed Directive Antennas 

m Determine beam shape , gain and number of beams 
to achieve those advantages apparent from sweeping 
antennas without the complexity of beam motion. 

Antenna Size 

• For those antenna gains found to be practical , what 
are the nominal sizes of antennas that provide! ihe 
desired gains and beam shapes? 
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3.1 


Sweeping Directive Antennas 


3.1.1 Sweep Rate , Gain and Duty Cycle 

Referring to the sketch shown In Figure 8, if a satellite is assumed to 
be at a radius R^ from the center of a spherical earth whose radius is R^ , then 
the satellite is in line -of -sight view of all platforms within a circle on the 
earth's surface defined by a cone with apex at the center of the earthy with 
semi-apex angle of c where cos c = R^/Rg • 

Having defined this viewing circle/ the total solid angle stibtended 
from the satellite to the circumference of this circle defines all directions 
from which platform transmissions can be received. If this solid angle at the 
satellite is called A^, then 

A^ = 2 ir (1-sln c) 

♦ 

This solid angle defines a reference point for comparing directive 
antennas of different gain. In particular/ an Ideal antenna with a circular beam 
whose dimensions Just cover the visibility circle would have a gain of 
As a numerical example, for a Nimbus satellite, this is an antenna with a gain 
of about 6 dB . 

For a single directive antenna of higher gain than 4Jr/A^ , complete 
coverage of the visibility circle will require a sweeping pattern. However, 
this pattern must be time and position ordered for random access systems . in 
particular, with platforms transmitting randomly In time from unknown positions 
within the viewing circle, the motion of the beam cannot be such that by chance, 
transmissions from a pai^lcular platform are missed throughout an overpass. 

More Importantly, where multiple transmissions separated by nominal intervals 
of time, are necessary, the sweeping pattern becomes further restricted/ordered . 
The manner in which these requirements can be satisfied Is by the Imposition of 
the following: 
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Satellite 
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• A directive beam must dwell in any given direction 
for a period of time greater than the interval between 
transmissions . This assures receipt of at least one 
transmission from each platform in the visibility 
circle during one complete sweep of the visibility 
circle by the beam. 

• The number of complete sweeps of the visibility circle 
must be no less than the minimum number of trans— 
missions required per overpass from the platforms near 
the horizon— i.e. , those with minimum viewtlme. 

Satisfaction of these two requirements Implies the following. The total view- 
time divided by the minimum number of transmissions <plus one) must be equal 
to the product of the number of discrete positions of the beam and the time 
interval between transmissions from the platforms. Analytically, these require 
ments may be written as: 

i = GA£T 

where 

Tv = period of time the platform is in view of the satellite 
during an overpass 

n = minimum number of transmissions which must be 

received from the platforms during the viewtime Ty 

G ® gain of the directive antenna 

~ (4ir/a) where “a" is the solid angle subtended by 
the beam of the directive antenna . 

T — ® Time between transmissions. 
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In writing these expressions, several simplifying assumptions are 
made. The ratio (A^/a) is the number of distinct positions of the beam neces- 
sary to completely cover the visibility circle. This means, there is no overlap 
of coverage from one beam position to another. Furthermore, the shape of the 
directive beam is such that each beam position is equally effective in cover- 
ing an area as any other beam position In terms of the projected area of 
platforms . 

For practical antennas, these assumptions are very optimistic in view 
of the relative geometry let alone the ability to shape beams arbitrarily. 
However, this assumption does place a lower bound on the number of beam 
positions and therefore total time required to completely sweep the visibility 
circle. 


By substituting, in the previous expression, the relationship for is 

namely 

At - 2.TrO-s«Me)- 2tr^\ - 1 - j 

The time between transmissions can be specified as a function of the 
gain of the directive antenna beam. That is: 



.^Tv 


(n+i) G 


- -^1 - ^ 


As mentioned previously, the minimum gain antenna that can be 
considered is that one which just covers the visibility circle which is 

r ^ 

Therefore, the time between platform transmissions that assures 
receipt of at least n transmissions during the viewtime T^ may be written as; 


_ Ty Gmin 
(n+1) G 
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Another bound on the gain, G, of the directive antenna can be derived 
from this expression. While is defined by the antenna that just covers 

the visibility circle, the maximum antenna gain possible is when the time 
between transmissions is equal to the duration of transmissions— i.e. , the 
platforms radiate continuously (CW). The possible spread of the antenna gains may 
then be defined as being between and the gain value corresponding, 

to CW platforms. A numerical example is interesting. 

For the Nimbus orbit, the minimum antenna gain is about 6 dB. If 
RAMS parameters are assumed, — l.e., Ty ** 5 minutes, n » 4, and transmission 
duration of 1 second— then the maximum (CW platforms) antenna gain is about 
15 dB. Therefore, the maximum advantage of directive antennas in the RAMS 
system for purposes of RFI suppression and/or decreased platform power is 
about 9 dB. 

3,1,2 Antenna Gain and Interference 

# 

While CW platforms establish an upper bound on antenna gain from 
the viewpoint of acquiring the requisite number of transmissions. Interference 
between transmissions occurring in the same beam may further limit antenna 
gain. To evaluate this interference, an assumption regarding distribution of . 
platforms is made. In particular, the platforms are assumed to be uniformly 
distributed as seen from the satellite. 


With this assumption, the number of platforms within the antenna beam 

of solid angle "a" is Na/A^ — where N is the total number of platforms within 

line-of-sight of the satellite. The rate at which transmissions from within a 

beam reach the satellite is then (Na/A„T), which upon substitution of the above 

A ) 

expressions for T, A^, and "a" leads to; 


Na 

TA„ 


N 


(n+1) 

— — = constant 

V 
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The conclusion, therefore, is that the time rate at which platform transmissions 
reach the satellite within a given beam is independent of the width of the beam 
itself. However, this can imply increased interference with higher gain beams 
as can be seen from the following. 


A rough approximation to the distribution of received frequencies at the 
satellite for platforms uniformly distributed on the earth is that these frequencies 
are uniformly distributed over a band, F,^.. If interference between simultaneously 
transmitting platforms occurs when they are closer than Af in frequency, '.and the 
transmission arrival times are pdisson distributed, then the probability of a 
given transmission being interfered with can be written as 


P = probability of interference = 1 - e 






where t Is the duration of platform transmissions. Nbte, the same assumption 


is made here, namely the number of platforms that can possibly communicate 
with the satellite is Na/A^— l.e. , only the fraction (a/A,^) of the total number 
of platforms as defined by the solid angle "a" of the antenna . 


At first glance, this expression would then indicate that the probability 
of Interference would be reduced in direct proportion to the gain of a directive 
antenna. This is true only if the other parameters remained fixed as antenna 
gain is varied. However, as shown in Section 3.1, the time between trans- 
missions must be inversly proportional to antenna gain in order to Insure receipt 
of a specified number of transmissions during a given minimum vlewtlme. By 
substituting this time and gain relationship, the probability of interference may 

be rewritten as: _ 4N(n+l)r Af 


Assuming the probability of frequency overlap, (Af/t^)) , remains the same, this 
expression then Indicates by the lack of terms describing directive antenna gain, 
that directional antennas will not Influence the probability of interference be- 
tween platforms * 
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Regarding frequency overlap. Figure 9 indicates that the probability 
of overlap depends upon the shape of the beam of a directive antenna as 
opposed to the gain. This figure shows, fora platform transmitting at 401 MHz, 
the locus of platform positions on the earth's surface giving rise to constant 
doppler shifts in received frequency at the satellite. Also noted are the areas 
encompassed by two different shaped antenna beams— one circular In shape 
and the other elliptical with its major axis parallel to the satellite subtrack. 

For the elliptical beam, the total range of doppler frequencies can be 
seen to be about + 7000 Hertz or just about the total spectrum of frequencies 
received if there were a single antenna giving complete coverage of the visibility 
circle. In this case then, the ratio Af/F,j,— the probability of frequency overlap- 
would not be expected to change significantly. Therefore, the probability of 
interference between platform transmissions would be independent of antenna 
gain as indicated above. 

This conclusion is not true, however, for the circular beam shown 
in Figure 9 . The total spectrum of received frequencies is on the order of 
+ 3500 Hertz or about one half the spectrum of frequencies for full coverage. 

In this case then, Af/F^ would be about twice the value for the elliptic or 
full coverage beams and the interference between transmissions would be 
correspondingly higher. 

From these considerations, the use of directive antennas implies the 
requirement to employ beams that are wide parallel to the satellite’s subtrack 
and narrow perpendicular to the subtrack. Otherwise, mutual interference 
between platforms will Increase monotonically with antenna gain. 

Because of the rather optimistic assumptions mentioned in the 
analysis of Section 3.1 , it is apparent that sweeping directive antennas 
cannot basically change mutual Interference for random access systems. 
However, other advantages are possible. In particular, the use of a directive 
antenna offers at least three opportunities to: 
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• Maintain platform power and transmission duration constant 
which, in effect, means utilizing the increased direc- 
tivity of the antenna strictly as a means to suppress 
interference from concentrated RFI sources— this might 

be advantageous even if losses associated with 
obtaining a directive beam negated any gain increase 

• Reduce platform power in inverse proportion to 
antenna gain in order to simplify the platform 

• Reduce the duration of transmissions from platforms 
in Inverse proportion to antenna gain by proportion- 
ately increasing platform data rates. This can be 
shown for random access systems of the RAMS type to 

t 

simplify signal processing on-board the satellite. 

3*2 Fixed Directive Antennas 

From the analysis of sweeping antennas, the possibility of reducing 
interference between platform transmissions is seen to be non-existent 
because of the necessity to decrease the interval between transmissions in 
direct proportion to increased gain. Furthermore, because of the optimistic 
nature of the assumptions regarding coverage of the visibility circle, a • ' 

sweeping antenna beam would inevitably result in performance degradation 
in that increased interference could not be avoided as the gain Increased. To 
alleviate this problem, an array of fixed antenna beams might be more advantageous 
However , the preferred orientation of these fixed beams may be dictated by con- 
siderations other than mutual interference. 

3.2.1 Beam Shape and Orientation 

The orientation and shape of fixed beams will be seen to be governed 

m 

by two considerations that are conflicting. 
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• From the viewpoint of sateilite signal processing, 
the antenna patterns should provide coverage such 
that transmissions are received in a narrow frequency 
band from each antenna 

• From the viewpoint of mutual interference between 
transmissions, the antenna patterns should provide 
full coverage of the received frequency spectrum. 

In essence, the shape and orientation of fixed beams is a problem of compro- 
mise between system capacity in terms of the number of platforms to be 
serviced and the complexity of satellite processing equipment— basically the 
complexity of the signal detection function on-board the satellite. 

To simplify the detection of transmissions, each fixed antenna should 
be elliptical in shape with major axis parallel to the zero doppler line— l.e. 
the major axis of the beam should be perpendicular to the direction of motion 
of the satellite. The minor axis of the beam’s cross section will then deter- 
mine the band of received frequencies from the platforms within the beam. 

From the relationship for mutual interference given in Section 3.1,2, this type 
of beam shape and orientation would not change the probability of interference 
between platforms compared to an antenna providing full coverage of the visi- 
bility circle— i.e.-, the antenna will accept all those transmissions that can 
potentially overlap in frequency. However, the detection processor’s lii^lified 
in this case because of the narrower spectrum which needs to be searched for 
arriving transmissions . 

On the other hand, to redude interference between Arriving transmissions 
the antenna patterns should reject in a geometric sense those platforms which 
give rise to similar frequencies at the satellite. This can be achieved by 
orienting elliptic beams with major axis parallel to satellite motion. If the 
separate beams stretch from one edge of the visibility circle to the other, then 
the exponent in the relationship for determining probability of interference will 





be reduced In near proportl^i^^i^^the. number of separate beams. However, the 
detection mechanism ass^diated with each beam must now cover the full spectrum 
of received frequenc^r^^ For example. If there are equivalently m comb filters 
required to detec|^nsmlsslons when a full coverage (single) antenna Is used, 
there will be cojnb filters required when there are b separate beams Involved 


V To analytically evaluate these two orientations of fixed beams, several 

assumptions analogous to those used to analyze the sweeping antennas will be 
made. First of all, the band of frequencies received from an antenna is 
assumed to be proportional to the angular width of the antenna in the direction 
parallel to the motion of the satellite. Secondly, the number of platforms 
within a given beam is, as before, assumed to be proportional to the solid 
angle of the antenna beam. 


With these assumptions, the probability of Interference may be 
approximated by 




\ - 


_ (5/1- e) 

TAt 


*v 


where 


Is the angular width of the beam cross~sectlon 
perpendicular to the motion of the satellite 


is the angular width of the beam cross-section 
parallel to the motion of the satellite 

so that A.J, is the solid angle subtended by the 

visibility circle of the satellite as before 
but approximated by a square in this 

^ . • C 

expression whose sides are (ff/2-c). 






is the fraction of the N platforms within 
the viewing -circle contained within the 
beam of the antenna 


is the fraction of the. total band of frequencies, 
F^, received at the satellite that will be con^- 
tained within the beam of the antenna. 


From this expression, the probability of interference is then only a function of 
the width of the individual beams in a dlrebtion perpendioular to the motion of 


the satellite—i.e. , the terms cancel. Therefore, for the case where beams 
have major axes parallel to satellite direction of motion and equal to (ir/2-e), 
the probability of interference may be written as . c 

I - e R- 


P 


where b is the number of beams. 


3*3 Size of Directive Antenna Arrays 

The half power beamwidth of a direct ^e antenna^ a^ay can be related to 
size by means of the approximate relationshiF< iJ= ('si AA)^ h^TirtiTe'l^lf'^ 
power beamwidth in degrees, Xis the wavelength of the operating frequency and 
L is the necessary antenna size perpendicular to the plane within which 0 is mea 
sured. With this relationship and assuming an operating frequency of 400 MHz ( 
meter wavelength) , the maximum antenna dimension is ,32b where b is the num- 
ber of antenna beams of the fixed type recommended in Section 3.2. Note , 
this last relationship is obtained by assuming a low altitude satellite (Nimbus 
or Tiros) wherein the total angle which must be covered by the b beams to com- 
pletely cover the visibility circle is about 120 degrees . As an example then, if 
there are three equal width antenna beams to each side of the satellite subtrack 

(i.e. , b = 6) , then the maximum dimension of the antenna arrays will be nearly 
two meters. 


Reference-Koelle-Handbook of Astronautlcal Engineering, pace 16-42 

McGraw Hill, 1961. y* y o 
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APPENDIX A 

PARTIAL DERIVATIVES OF THE MEASUREMENTS 


Section 2.4 shows how a platform may be located using an iterative 
technique for solving certain non-linear equations. Section 2,5 showed how 
the mean location error may be computed from the measurement errors , Both 
of these computations Involve the use of the partial derivatives of the mea- 
surements with respect to the trajectory parameters (e.g., equation (5) of 
Section 2,4 and equation (1) of Section 2.5) . These partial derlvates are 
computed in this Appendix. 

A.l With Respect to Transmission Instabilities 
A . 1 . 1 Doppler 


In a Doppler location system, the received frequency (fj^) is measured. 
The relationship between the received frequency at time tj^ and the trans- 
mission instabilities is 



where fti is the frequency transmitted at time t^, c is the speed of light, and 
is the time rate of change of the distance between the platform and the 
satellite at time t^. The transmitted frequency, f^^, is approximated by 

where ^to is the nominal transmission center frequency, f^g is the center fre- 
quency drift, ft is the centeh frequency drift rate, and t© is a reference time. 
The reference time Iq may be chosen as the time of the first transmission on 
an overpass. In the case of a multiple overpass geometry, the frequency drift 
and drift rate may then be computed for each overpass. 


The partial derivatives of f|^^ with respect to f^g and f^ are 





■ I - 




c. 




) 




A. 1.2 Doppler Rate 

In a Doppler rate location system, the time rate of change of the 
received frequency (f|^) is measured. The relationship between fj^ at time 
and the transmission instabilities is 


V " C / y 7 


a 


dVri 


where is the radial acceleration between the satellite and the 

platform . Then 


P/iT 




^ o 


<±n - ( ^ ^ I 

gp t C J ^ ^ ^ C. 


A. 1.3 Range Difference 


In a range difference location system, the difference in time (At) 
between received transmissions is measured. The relationship between At^^ 
at time ti and the transmission instabilities is 




^^/(c ^C. 


/?.3 


where ARj^^ is the difference in range at times t^^ and and Nq is the number 
of cycles between transmissions. Then 
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9 





K/c 



9 h 



(U - i,) 



A. 2,0 With Respect to Platform Kinematics 
Define the following symbols; 

1) Rg^, Vg^, Ag^ - platform position, velocity, and acceleration 

vectors at time t^ 

•—A » ^ ^ 

2) Rgj^, Vgj^, Ag^ - satellite position, velocity and acceleration 

vectors at time t^ 

3) Rj^^, = relative range, relative range rate, and 

relative range acceleration at time t^. 


Platform position, velocity and acceleration at time t^ are functions 
of the Initial state of the platform, i.e. , its position/velocity/ acceleration 
at time 

~ ^ r ^6® i ^ Oo ) 

1 

^et ^ i' i ^ao j ^ e<, j l 

Each component of all vectors functions 

of each component of the vectors Rg^, Vg^, and Ag^. It is the initial state of 
the platform that is to be estimated. It will be convenient to use the notation 
* to denote a component of the initial state of the platform. For example, 


would refer to the partial derivitive of platform position with respect to RBxq 
when * = , with respect to Rb^ when * = Rgy , etc. 

A. 2,1 Doppler 


The relative range rate is 


V = ( ~ ) 




The partial derivatives of are 


3 Vf: 

3 * 






k e r c 





'~rUcn 




£lL 

o ^ ^ 
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A . 2 .2 Doppler Rate 


The relative range acceleration is 


^ 1 . 








The partial derivatives of are 


S>^ 


I . 




9^ 


(^si" 6l) 






then 


# • # 


9/e 

^ ~ £jl - 


9 fit? ^ 



C 

9>if 

- 

.0 

'2'. - £l P ^f^t 


9 fie: 


c- 

d. 



^ ^ 1 

A. 2, 3 Range Difference 

The relative range at time t^ is 

^OC I 
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The range difference is then 


^ =r / ~ 


The partial derivatives of are 


3 


9 ^ 0 / 


=> — 


then 


5 ^ 





APPENDIX B 

MODEL OF PLATFORM MOTION 


This Appendix presents the model of platform motion used in the 
TRILOC location error analysis . In this model of platform motion, the plat- 
form is at constant altitude with an initial horizontal velocity (which may be 
zero) and an Initial horizontal linear acceleration (which may also be zero). 

This model approximates a great circle trajectory. The following symbols 
are defined as: 

V, A - platform position, velocity, and acceleration 
at time t. 

^o'^o “ platform position, velocity, and acceleration 
at time t^. 

6t - t-to 

C - platform altitude (constant) 

Then, 


- 


(Si)^ 



- 

Si ^ 

r it) ^ 




i C-- V- f - 




V/j 



^*3 ~ 

0 




/?e 


- 

^ ^ V- n 


/9y- = - 

VX ¥ ■/ ■* Pi 




' X 

Px 


B-1 



^ The partial derivatives of P, V, and A with respect to the components 
of Pq, Vq, and Aq are 


9 P -■ 


3 P 


y e> 


- ::6 - i : 
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.5 P 
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^ ^ ^ i 
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_ ^ ^ X'1^ ^ 

= 
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9P 
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*7'^^=^: -:i f rr^7' '?'o 


P>f 
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<s^ 


-Y 


i-..- ..- pp 




p ^ P/y 

^ t ^ ^2 ^ 


^ . (f" ^X. 


The derivatives of P, V, and A with respect to P^o# V^o-, and Axo at® zero. 
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APPENDIX C 


n DIMENSIONAL HYPERSPHERICAL 
COORDINATE TRANSFORMATION 


In Section 2.5 the expected location error was written as a multiple 
integral. It was stated that by performing a transformation from Cartesian 
coordinates to hypers pherica I coordinates the integration could be performed 
over all variables except one, i.’e. , the multiple integral could be reduced 
to a single integral. This Appendix defines the hypers pherica 1 coordinate 
system and evaluates the determinant of the Jacobian of the transformation 
from hyper spherical coordinates to Cartesian coordinates. 

Define a transformation from n dimensional hyperspherical coordi- 
nates (r, 9j, . . 6n-i) to n dimensional cartesian coordinates (x^, . . .,Xp) 

as follows: 


n -/ 


Xt ^ 7T ~ . . C OS. _ 

1 . 


n 


77 ". c^s6^/c . _ 

c 


L ^ A. A ^ 




c-i 


The inverse transformation is 



The range of the hyperspherical coordinates is 

■ o 


-. 77 !?. srr 

.'r.jry X. -.ir/'t- .. 

The Jacobian Jj^ of the transformation is the matrix 


^ C ^ * * "i X’rt ^ .... 

" V 



cP^, 


. , 

♦ , , _ ^ 


w 

• 

V 

• 

» 

rr 
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We want to show that the determant of Jn is 





I- - ^ , /7 

iss/ / 


For n = 3 , we have 


X, 

— X’z -- irt . tf'j i 

A'3... *_..r s tin, . ... 

and 

_ 6ff CoS 6^z T'X.stn ^f XoS.6^Xr- ^ tCo: <>, sT//7 '<?^ 

A^J / .. . .Hrt.A/. CoS Sz tc ..CoSt^fCOS^z 

... j... s±r> <J? CoS &Z, 


r . c.a^ &z.. 


r ir^ 




C oc 


y-y-- 


Therefore, the formula holds. for n = 3. 
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Now suppose 


I 'Jfl./ I ~ 




0-3 


" TT- i 

r JT cos 

I't 


for n-i 3 , Define 


/I ^ ^ 

f ^ . C 7T 

/T* z. 

r . rs/r) TT, 

A*^ 4 




~s— Then x^ = yi cos 9n-l for i ^ n-i and jCn = r sin 6jn-l . The transfor- 
mation from (r, 9i, . . ., 8n-2) coordinates to y^-l^ coordinates 

is an (n-1) dimensional hyperspherical coordinate transformation . Therefore, 


' S>( ?/j 


^ C Tj ^ > i 






~7T^ c &j 
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If we let zr 



then for i ^ n-1 , 


r- 


^ '( iJc' ~V>) <;> t* 


^ r 







Cos 


~ ~B Xc ^ ( 'i c COS^n--'^ ^ 

S>6>J-r -- --^6^r- 


Cos 


- -. =! - ^ I S- - CoS 6'r)-t ... C .?. . ^ . j' 

* • ' 1 

Qx: _ ^ i'i: Cos 


cP 6fn^', 




- 'T..IC...^ i i ..-5//?, _ 





and 

« 

^ /- 








c? / ^ S ^ r) — 

-X CoS 4^/7-, ^7" J ~ /7 
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We then have 


a.,/ = 


^ r7.f. / 6*r7-/. 

t C«S&/f n~i_Coi. 


/»- / - /. ._ 

. ..•£>>? 1 ■ ■ * ■* - * — ^ — - 
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,__=Y-/; 




^ ^ CoS — ! ^/ 7 -/ Cos On^ i 
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# ■ 

^(V-/ . t ^/7~0 ^ n-.^'L^CiS^rt'^i r 0 * 'c&s 6^r>'‘ 

' ■ ''- ^ _„<?. - „ .(? 


of the 


_ ^ /»? ^/? -/ rcos On - 

\ 

This determinant is easily evaluated by observing that the cofactors of 

/i 1 ^ no ■rt ^ t<r +-]o^i 1 o 


first two elements of the last row are 
1) the cofactor C\ of sin 9^-1 is 


! 

“ ^ ^ i X - ^CS, I .. _ • . /7^ t.. CCC 

V (s"0 I S/^^02^/ . ^. , CM 0 


^ .^/i i Sin 6 ^n 


C^S 4 C^S Cn^ j 


C-6 






S Co 5 


C,J 


2) The cofactor of r cos 0n-l is 

t 

: . ^/X . CoS — Z 


^7* Z 


i' 






I Cfijj£>n-!^ . C9S^„^, ., • ... 



r7yi 


“ (^0 ' / '^/»-/ / /* ^ 


The determant of Jj^ is 

iJnl *“ ^/ / “/■ ^ i. 
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We then have 


‘ - /7--2 


’ ~ I fi-* Z 

IT, CoS C'~0 L'^n^i- 1 — ^/7 jl- 


1 7" ’ X 
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^ ^ J -V 
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^ 7^ - 
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jr._A.<^A..^ Uj. 


Thus w^ have shown by induction that 

IJrt / ~ r'^ * It ecu 
Isr-3 

for all n < 3 , 




APPENDIX D 


7^/t 


Evaluation of 


/ 

-TTA 


F-f , 

Cos 




2./S s,r^o c«xe^ V- y 


t 

Sf*^ 


X *7 


The mean location error was written as a multiple integral in Section 
2tS. In Appendix E the multiple Integral Is reduced to a single integral. 
Beginning with the Integration over each leaves an integral of the 

form 


where a, and X are functions only of , and 6, The Integral D-1 . 

J-z 

is evaluated in this Appendix. This result is then used in Appendix E to com- 
pute the mean location ezror. 


Note that 

7r/^ • _ 

Cfii e/<p 

^ ^ / 

J 6^ ~-jh ^ V y “ 

7T/^ ^ 

CgJ ' 

o 

/•-/ ’ 

r Ce>i fg* 

'^ J [o(r^ctfS^6>’~- 2./S s/r>o y 

For convenience, we shall denote the integral to be evaluated by i (e). 
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Then 
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M “ and v 


when t = 0, and u = v -o^ 
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APPENDIX E 


Evaluation of 


Tr/"*- 


TT/l. 






-77/2, -TT/2 


r . . , 


In Section 2,5 the mean location error was written as a multiple 
integral. This Appendix shows how the multiple integral is reduced to a 
single Integral. No analytic expression for this single Integral was found 
but its value can be numerically computed. 


. P = _ . fiT.tlA.. . - 




= r Cc^S^j *77“ 1 ♦ i * '-V ^ 

-L-- • ^ 1 rt»t — 


& 


o^t 


Thus P is an n X n dimensional matrix and y is an n dimensional column vector. 
Define 


and 




E-l 



Then for 


^ ^ Z, /} f 




and 


2 . ^ — 4 * 7 ^ ^ 

p;., ,.= Vj*. Pj-..~^ rj r/-. 




where 


Pi- 


p- .. = _^_. 


' r, 


I y; 


★ 

That is, Pj is a j X j xliinensional matrix. Pj is the (j-1) x (j-1) dimensional 

matrix derived from Pj by striking out the row and the column of Pj . 

?j is a (j-1) dimensional row vector which consists of the j row of P^, but 

without the (last) element of that row. '/.• is the jfl' element of the 
“* * 

row of ?y' Yj is a j dimensional row vector such that 

1) ?n - y • 

2) yj«i is computed from yj by striking out the last element 
of yj and dividing the remaining j-1 elements by cos . 

Dashed lines in the expressions above are- used to indicate partitioned matrices 
and partitioned vectors. 
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Then, if 3 s j in, we have 



= ../... . 5 ;-- P J 'iS-> J --' 



*[yj]. 


^ z 

,_"z J (?j 


where i 







Note that*:^® and ^ ^e scalars. 
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We then have 



We now want to show that 


— ir/x. TT/v,, r.j!2z' r n . , 
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Let n = 3 . Then (see Appendix D) 
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Therefore , the formula holds for n = 3 . Now suppose it holds for n - 1 a 3 . 
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“ . 2 p(-^) [ j- 

Thus we have shown by induction that the.fonnula holds for all n ^ 3 , 
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APPENDIX F. 

MEASURING FREQUENCY AND FREQUENCY RATE 

This appendix describes one technique whereby the frequency received 
by a satellite and its time rate of change can be estimated and determines the 
precision to which these two parameters can be measured. The basic assumption 
throughout the analysis is that white noise is the only error source corrupting 
the measurements. 

The measurements from which frequency and frequency rate are esti- 
mated are the series of elapsed times between positive going zero crossings of 
the received frequency {after down converting the r.f . signal to several thou- • 

sands of Hertz) . For present purposes, the resolution and precision of these • 

# 

time interval measurements is assumed to be much smaller than the equivalent 
phase jitter of the signal itself as caused by noise. 

Assuming the received frequency is varying linearly with time, then 
the period of time (T^) between i positive zero crossings can be related to 
frequency (w) at a given point in time and the time rate of change of frequency 
<a) by 


F-1 



If perturbations are taken from this equation, then the errors in a, w, and Ti 
(6a, 6w, 6Ti) are related by 

^ Sot — C'^ 

or for an ensemble of measurements of T, this may be written in 
matrix notation as 

where for the ith measurernent 



With k>2, the errors in 6w and 6a can be determined if the solution 
for w and a is assumed to be based upon least square fitting or regression. 

In this case, the solution is 

(a"^a)‘Vst 

By assuming the 6 t's are independent (which they are not strictly speaking) of 
zero mean and equal standard deviation of ff-p, then the covariance matrix of 
6x becomes 

cov (Swj?pc) ^ 
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This relationship can be used then to evaluate the standard deviations 
of w and x (namely 0^ and (To,) as a function of w and oc. This has been done by 
means of a computer program with the following results. 

Wiih frequencies (w) and frequency rates (oc) between 1 to 10 kHz and 
2 to 100 Hz/sec, the standard deviations of frequency and frequency rate are 
essentially independent of frequency rate and approximately proportional to 
reciprocal of the square root of the frequency itself when measurement duration 
Is fixed. Fora duration of .5 seconds, 

0^-Js/k^ ^ .07 V\% 

- .% 7 /^ 

where f is the nominal or average frequency in units of kilohertz and cr-j. is 
replaced by the relationship between signal to noise ratio and frequency — 



for w in radians/sec 
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